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We study a simple quantum mechanical symmetric donor-acceptor model for electron transfer 
(ET) with coupling to internal deformations. The model contains several basic properties found in 
biological ET in enzymes and photosynthetic centers; it produces tunnelling with hysteresis thus 
providing a simple explanation for the slowness of the reversed rate and the near 100% efficiency of 
ET in many biological systems. The model also provides a conceptual framework for the development 
of molecular electronics memory elements based on electrostatic architectures. 
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The physics of electron transfer (ET) is of paramount 
importance in many fundamental aspects of molecular 
biology, from photosynthesis to mitochondrial ET, to 
many control reactions and signalling through the cell 
membrane |3, |3, Being a sophisticated quantum me- 
chanical tunnelling process in itself, the ET problem in- 
cludes many aspects of coupling to vibrational degrees of 
freedom, tunnelling control through external variables, 
and non-reversibility aspects ,2]. After the pioneering 
work of Marcus p| in the electrochemical aspects of ET 
for redox reactions in solution, the problem has seen 
several times renewed interest in both semiclassical and 
quantum mechanical (QM) aspects. This is an area 
where there is a strong interplay between QM and im- 
portant biological effects at the molecular level. 

In addition, there is considerable evidence now on the 
importance of nonlinear contributions to many impor- 
tant dynamical aspects of molecular biology. DNA dy- 
namics and denaturation[(| Q and signal transduction 
and coordination of events during the DNA readout by 
proteins^, [9| have both been ascribed to nonlinearities. 
The nonlinear character of the hydrogen bonding poten- 
tial a mong DN A bases has been known for a very long 
time pi flOl fll| ■ But nonlinearities in biology are believed 
to play also a substantial role in ET and energy transfer. 
Resonant dipolc-dipolc interactions of the amide I band 
combined with structural deformations in alpha-helical 
proteins have been proposed for a very long time as a pos- 
sible explanation for energy transport from the hydrolysis 
of ATP [!2 • There is experimental evidence that amide 
I vibrations are coupled nonlinear ly; an example being 
their unusually long life in simple proteins|l3j. The cou- 
pling of carriers to conformational degrees of freedom has 
been suggested also to exist in connection with ET from 
the active prosthetic groups inside cytochromes Q. There 
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seems to be, in addition, substantial evidence that a non- 
linear coupling is required to explain the so-called gated 
ET reactions between quinones in the photosynthetic re- 
action centers [T^. Knapp and coworkers |l5j were among 
the first to realize the importance of both nonlinearities 
and Davidov-like solitons 4] in the ET problem. They 
argued that the construction of the initial and final wave- 
functions in terms of solitary states leads automatically 
to localization for a symmetric donor-acceptor system, a 
very desirable model property in many cases, as we shall 
explain later. 

The brief outline given above highlights the fact that 
there is a case in many real biological systems for the 
study of QM tunnelling properties in the presence of in- 
trinsic nonlinearities. To this end, there are two possible 
options, either we study the specific non-linear coupling 
of a given case with all its peculiarities and details, as 
done in Ref. 0] for ET between quinones in the pho- 
tosynthetic reaction center, or we study a skeletal model 
which contains the essential elements of the problem and 
try to draw general qualitative conclusions from its physi- 
cal behavior. We shall adopt the latter view here to high- 
light a specific property of the tunnelling process in the 
presence of nonlinearities; i.e. the presence of hysteresis 
in a symmetric donor-acceptor ET system. 

Davydov's approach to electron-lattice deformation 
coupling in biological systems^] and its variants^ lead 
always to an effective Hamiltonian of the non-linear 
Schrodinger equation (NLS) type. The nonlincarity in 
the equation comes from a hidden adiabatic coupling to 
an external interaction (the deformation) which is fac- 
tored out from the dynamics in the form of an effective 
self-consistent potential. The Hamiltonian is always of 
the self-localizing type, i.e. with a negative nonlinear po- 
tential. We have study very recently two models with this 
property in the framework of the ET problem [16] . The 
simplest possible discrete symmetric donor-acceptor ET 
system, including a nonlinear coupling of the Davydov 
type, is the two-sites nonlinear Hubbard (NLH) model, 
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given in second quantization by the Hamiltonian: 



H = a\c\c\ + e 2 c\c2 + t(c\ci + 
U (ni) cjci + {/ (n 2 ) 4c 2 
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where (rij) = ^cjcjy, and 1 (2) refer to the donor (accep- 
tor). The Hubbard-like terms oc U > give the neces- 
sary (non-linear) self-trapping potential and Ci and t are 
the site energies and overlap (hoping) integral between 
the two sites, respectively. This would be the case of 
tunnelling between two equivalent sites with strong po- 
laronic coupling. This model, and variants of it including 
dissipation, has been extensively studied for many years 
by Kenkre and cowor kers[3 13 13 13 who coined the 
term quantum nonlinear dimer. The objective of this 
Letter is to focus on a specific property of the tunnelling 
between the two sites in this model, namely: the presence 
of hysteresis. We shall then speculate on situations where 
this property may be playing a role in ET with nonlinear 
couplings. We review very briefly a few of its properties 
discussed in Ref. 0. The ground state can be sought 
by a Monte Carlo (MC) simulated annealing in \z\ — z 2 ] 
space, were z% and z 2 are the complex coefficients of the 
ground state wavefunction for the amplitude of the elec- 
tron at the donor and acceptor sites, respectively, i.e. 
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with \z\\ -\-\z2\ = 1. A good measure of the asym- 
metry in the electronic density of the two sites is the 
magnitude 



S = 



P2 
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which gives the absolute value of the difference between 
the probability densities for donor and acceptor. Figure 
^is the calculated S as a function of t/U in the ground 
state of a symmetric donor-acceptor system with e± — 
(-2 = 0. Alternatively, the minimization of Q subject 
to the normalization condition renders S = Wl — (t/U) 2 
for t/U < 1 and S = for t/U > 1, as seen in Fig. [I] If 
t/U is large (small U) the ground state is symmetric, i.e. 
S = 0, and the electron is equally shared between both 
sites. Below a certain value of t/U, the system displays 
a spontaneous symmetry breaking in the ground state, as 
discussed in Ref. From there on the system profits 
from accumulating the wavefunction in one of the sites 
and S — ► 1 as t/U — > 0. This is the localization in the 
symmetric donor-acceptor system mentioned by Knapp 
et. al.[3- 

The question we raise now is how the electron behaves 
when we start varying the energy e, of the levels exter- 
nally. We shall show that the intrinsic nonlinearity of 



Davydov-type models for ET introduces hysteresis in the 
quantum mechanical tunnelling. 

The qualitative reason is relatively easy to understand 
and is shown in Fig. [21 Let us assume that U is big (t/U 
small) and the electron is mainly trapped on the left. We 
now vary the potential of e 2 making it lower than e\. 
Under normal circumstances, the electron would tunnel 
immediately to the right. However this might not be the 
case in the presence of a nonlinearity, because the energy 
of the electron localized on the left is of the order of ~ 
e\ — U, and this could be less than e 2 , i.e. the nonlinearity 
is holding the electron in its original place. There is a 
memory effect, accordingly, which cannot exist in any 
linear, adiabatic, quantum mechanical system. The same 
holds for the reverse situation if we raise e 2 above e\ once 
the electron is localized on the right. By fixing ei, the 
tunnelling events left-bright or right— deft occur at two 
different values of e 2 , producing an intrinsic hysteresis in 
a very simple QM system. 

Mathematically, however, the property of hysteresis 
with a nonlinear Hamiltonian is in general difficult to 
prove. One could argue that in Fig. |3 a gradual spread 
of the wavefunction from left to right is produced un- 
til the electron is fully localized on the right. And the 
opposite might hold on the reversed bias situation, thus 
producing a crossover without hysteresis. The way to 
prove the existence of hysteresis with a sudden transi- 
tion between the two situations is by properly solving 
the model. Notwithstanding, a perturbative approach is 
doomed to fail. The reason is that a localized transition 
from donor to acceptor is more of a cooperative effect. 
In what follows we concentrate on the problem of the 
stability of the wavefunction given an initial (localized) 
starting condition and a bias between t\ and e 2 . 

It is advantageous to change to a continuous version 
of the model in Fig. |2 We treat a model with a sim- 
ilar property (spontaneous symmetry breaking) that is 
related to the nonlinear particle in a box. In Ref. 
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it was shown that a quantum mechanical particle in 
a box obeying the NLS equation undergoes an sponta- 
neous symmetry breaking in the probability density for 
a sufficiently large value of the nonlinear potential. The 
particle is confined either on the left or the right of the 
box and it is possible to produce tunnelling between the 
two situations under the presence of an external stochas- 
tic potential. The NLS, which is the equivalent to for 
a continuous system, reads 



7 VV(r)+/3|^(r)| 2 V(r) = i?V(r), 
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where ip(r) is subject to the boundary conditions 
"0(l)l±a/2 = 0, with a = size of the box. E is the energy 
and 7 < 0, [3 < are the kinetic energy and self-trapping 
potential parameters, respectively. The relevant magni- 
tude is the ratio p = \j/0\, fixing the relative strength 
of both. For small p, an spontaneous self-localization of 
the wavefunction either to the left or the right of the box 
occurs [l^, both sides playing the equivalent roles of the 
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donor-acceptor sites in . The continuous model J3J al- 
lows the inclusion of an electric field competing with the 
stability of the self-localized ground state. The presence 
of an external potential, in the one-dimensional problem, 
is added to J3J through an additional term on the right- 
hand side of the form: 



(5) 



where Vq is a constant, and r = x in this case. This 
potential produces a linear slope (constant electric field) 
form one end (x — —a/2) to the other (x = +a/2) of the 
box. 

In addition, the stability of the wavefunction can be 
tested with the aid of the Fourier Monte-Carlo (FMC) 
method developed in Ref. 0. We review here very briefly 
the basic concepts involved, while leaving details to the 
references ^(J. The wavefunction can always be expressed 
by a Fourier series of arbitrary length of the form 
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with k* = 2n/a (n + 1) and k^ = 2n/a (n +1/2). a n -b n 
are Fourier coefficients and ^(r) satisfies automatically 
both the boundary and normalization conditions. The 
basic idea of the FMC method is that the ground state 
is found through a simulated annealing (as a function of 
an artificial temperature T) in Fourier coefficient space 

Q"n~bn ■ 

Once the ground state is found for T — 0, the poten- 
tial can be changed and the evolution of the ground state 
monitored for small finite temperature To- It is under- 
stood that the particle follows adiabatically the changes 
in potential. The role of T is that of a stochastic per- 
turbation testing the stability of the ground state; T 
has to be smaller than the energy of any metastability of 
the system in the presence of nonlinearities. The conver- 
gence of the method is tested by doubling the number of 
coefficients and comparing the quantitative results[lfij. 

Figure 13 shows the result for a full cycle of the electric 
field in both directions. The parameters of the simulation 
are given in caption. We start with the particle confined 
on the left (which we arbitrarily call donor) and we follow 
the evolution of the wavepacket by monitoring (x) /a = 

a/2 

(I /a) J xtp 2 (x)dx for different potentials in the cycle 

-a/2 

shown in Fig. |^a) . Figure D ) shows a result with 
the hallmark of hysteresis. As before, the reason is the 
competition between the intrinsic nonlinearity and the 
external potential, fixing either the mismatch of the levels 
in Fig. [5] or the energy difference between left and right 
in Fig. IHIb). The transition is not gradual (as it would 
be without nonlinearity) because the electron does not 
profit from a reduction in energy by transferring part 



of its wavefunction to the other side, until a very specific 
condition in terms of the competition with the correlation 
energy U is achieved. 

By having a model that combines basic aspect which 
are know to exist in real ET systems, we can now go one 
step backwards and contemplate those situations where 
this effect might be playing a role. Control of ET is 
decisive in many biological functions; the canonical ex- 
ample being ET in the mitochondrial chain, which must 
stop when ADP is all phosphorylated. For more than 
30 years it has been suggested that the energy matching 
between donor and acceptor could provide a mechanism 
of ET control[^l In general terms, this ability of 
control has to be compatible with the slowness of the 
reverse rate responsible for the near 100% quantum effi- 
ciency observed, for example, in the photosynthctic reac- 
tion centers, and with the presence of couplings to the lo- 
cal environments. These requirements are normally met 
by the so-called energy-gap law which postulates two dif- 
ferent activation barriers for the forward and backward 
ET reactions, in the spirit of Marcus theory 0, 0, HJ. It 
has been also known for a long time that simplified in- 
vitro model systems 23J invariably show similar or larger 
rates for the backwards recombination. It seems then 
that vibronic nonlinear couplings to the environment are 
essential to understand the efficiency of the ET mecha- 
nism, and that evolution through natural selection has 
engineered a favorable combination of nonlinearity and 
control capabilities. 

It is clear, accordingly, that a basic model like the one 
presented here satisfies several of these conditions simul- 
taneously. The presence of QM-hysteresis is a key issue 
that has been overlooked in the opinion of the present 
author, and provides a very fine mechanism of control 
even in the presence of broadening of the electronic states 
through vibronic coupling at the donor and acceptor sites. 
In addition, the back-reaction problem is automatically 
solved by the same token, without having to resort to 
ideas which are normally more difficult to justify, like 
the presence of two radically different matrix elements for 
the forward and backward reaction, respectively |24(. We 
believe that a model of tunnelling between donor and ac- 
ceptor in the spirit of the models proposed by DavydovQ 
and Knapp[15|. contain the essential qualitative ingredi- 
ents to explain simultaneously both control and efficiency 
of real ET system found in biological systems. The main 
objective of the present paper is then achieved in Fig. 
13 i.e. to demonstrate the existence of QM-hysteresis in 
a very simple model and to highlight its possible rele- 
vance in a variety of phenomena associated with ET in 
biological systems. 

The hysteresis cycle in Fig. [^defines a threshold value 
for the transition between the two states, in analogy to 
the coercive field displayed by magnetic systems. In or- 
der to show that the hysteresis in the cycle shown in Fig. 
13 comes effectively from a competition between the in- 
trinsic nonlinearity of the Hamiltonian and the external 
electric field, we show in Fig. 21 the effect of the param- 
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eters of the model on the coercive filed of the transition. 
Figure 0] shows three examples for a fixed 7 and three 
different /3's: -0.1, -0.05, and 0. For = -0.05 the cy- 
cle has narrowed down to a very small region around the 
origin. The system still profits from accumulating the 
probability density on the left or the right of the box, but 
jumps immediately form one situation to the other un- 
der the slightest perturbation. The Fourier-MC method 
has problems to define the exact width of the transition 
in this case, because the perturbation parameter T can 
be comparable to the small barrier separating the two 
situations. For f3 = 0, the particle follows smoothly the 
external electric field. The S-shaped average coordinate 
of the probability distribution is due now to a competi- 
tion between the boundary conditions on both sides and 
the force exerted by the electric field. The wave-packet 
follows the field linearly at very low fields, but as soon 
as the displacement from the center is large enough, the 
boundary conditions on the left or right of the box limit 
the response to the field; thus achieving saturation in 
the displacement. These results prove that the nonlinear 
bistable behavior in Fig. [3] is, effectively, a competition 
between nonlinearity and external fields in the model. 

On more speculative grounds, however, the existence 
of simple QM-systems with tunnelling hysteresis at the 
molecular level should also be of interest in other prob- 
lems. For example, in the field of molecular electronics 
(ME) [2^] where decisive steps have been taken towards 
the production of logic devices based on the change of 



electronic properties by single molecules. Charge storage 
in a self-assembled molecular device has been demon- 
strated already!^; the device works by storing a high 
or low conductivity state achieved by the addition or 
subtraction of an electron. But the contact problem 
of organic molecules with metals remains an outstand- 
ing problem in ME, and alternatives have been proposed 
based on the so-called quantum cellular automata (QCA) 
and electrostatic architectures |25l |27| . The QCA ap- 
proach has received considerable attention in ME as an 
alternative to conventional electronic devices due to its 
potential in reducing power dissipation by avoiding exter- 
nal currents [2j| . It is also an interesting conceptual shift 
from conventional electronic design. After all, the most 
sophisticated machines known in nature, which are living 
organism, work mainly not through changes in conduc- 
tivity and currents but rather through electrostatic and 
van der Waals interactions. 

Even though the technology to achieve ME and QCA 
devices is only beginning to appear, the results in this pa- 
per suggest that organic symmetric donor-acceptor pairs 
with strong polaronic coupling, and external control by 
an electric field, are excellent candidates for the ultimate 
single-molecule dynamic random access memory (RAM) 
device. 
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FIG. 1: Asymmetry parameter S (Eqn. [3J for the ground 
state as a function of t/U for the two-site NLH model (Eqn. 
0. A schematic view of the model and its parameters is 
given in the inset. The calculation is for a symmetric donor- 
acceptor pair with t\ — £2 = 0. Below t/U = 1 the system 
profits by accumulating the electron on one of the sites, thus 
creating an asymmetric self-localized ground state. 



FIG. 2: By assuming that the electron is initially localized on 
the left (as shown schematically by the plot of ip 2 on the top), 
tunnelling to the right is hindered even if 62 < ei because of 
the correlation energy U. The system has a memory of its 
ground state producing a hysteresis loop for the tunnelling 
rate between both sites. See the text for further details. 



-4|- 

0.0 5.0x10 5 1.0x10 6 

MC-steps (time) 




-4 -2 6 2 4 
V [Volts] 



FIG. 3: Hysteresis loop for the particle in a box with a nonlin- 
ear potential (Eqn. In (a) the amplitude of the potential 
is shown as a function of MC-steps, which play the role of 
time. The parameters used in the simulation were 7 = —1.0, 
(3 = —0.1, and the size of the box a — 0.5. Further details 
of the numerical implementation of the FMC-method can be 
found in Ref. ITrl The particle follows adiabatically the po- 
tential changes along the MC-changes. The initial state is 
with the particle on the left, and it is found after a simulated 
annealing with 20 temperatures from T = 1 to 0, with 10 4 
thermalization cycles per temperature and 40 Fourier coef- 
ficients in total p^. To is fixed subsequently to 0.05 during 
the potential changes. In (b) the evolution of the expectation 
value of the wavepacket {x) /a is plotted as a function of the 
external potential. Note the presence of two very different 
thresholds for the tunnelling in both directions. As in Fig. 
the reason is the competition between the intrinsic nonlinear- 
ity of the Hamiltonian (once the particle is localized on the 
left or right) and the applied external potential. 
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FIG. 4: Hysteresis loops for different values of the nonlinear 
coefficient j3 in Eq. 21 For a fixed 7 = — 1.0, we show three 
examples of competition between nonlinearity and the exter- 
nal electric field for j3 = —0.1 (top), = —0.0.5 (middle), and 
/3 = 0.0 (bottom). For /3 = —0.05 (middle) the coercive field 
on the top graph has collapsed to a very small value around 
V ~ 0. The system profits from accumulating the probability 
density on the left or the right of the box, but the barrier is 
very small and quickly overcome when the field reverses. For 
P — 0.0 the particle follows (linearly) the electric field at low 
voltages, until the boundary conditions on both sides of the 
box start competing (and saturating) the drive of the external 
field. See the text for further details. 



